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Abstract. Arguments for H and B to be considered 'axial', or pseudo vectors, 
are revisited. As a point against, we examine the complex-coordinate method 
for numerical grid truncation and mode loss analysis proved very successful in 
computational electrodynamics. This method is not compatible with convention 
that H and B are axial. 



1. Introduction 

In physics, one often encounters 'symmetries', i.e. situations when structural fields (as 
the permittivity e and permeability /i in classical macroscopic electrodynamics) are 
invariant under the given transformation x' = x' (x) . To obtain general conclusions on 
dynamic behavior of such symmetric systems without solving the governing dynamic 
equations, one needs to know also how the functional fields (electric and magnetic 
fields in Maxwell's equations) are transformed with x' = x'(x). Physical intuition and 
experiment is generally what provides us with this knowledge. 

Important symmetries are with respect to coordinate reflections, as x — > x' = 
(— x , x 2 , x 3 ), and coordinate inversion, x — > x' = (—a; 1 , — x 2 , —x 3 ) — — x. These are 
improper transformations whose Jacobian, J? = dx" 1 jdx 1 ', has negative determinant: 
det J\ = A < 0. In numerous textbooks on classical mechanics and electrodynamics, 
vector-like objects are classified into polar, or ordinary vectors, and axial, or pseudo 
vectors — depending on their behavior under improper transformations. Polar vectors 
are transformed under inversion as V = —V, while axial as W = W. 

The V = —V transformation under coordinate inversion implies that the V 
vector remains unchanged in the reference ('absolute', 'physical') space. On the 
opposite, W is reflected in the absolute space upon coordinate inversion. How is 
it possible at all that such 'unrealistic' pseudo quantities, dependent on the way we 
describe them, survive in the physical picture of reality? The answer is that they are 
normally related to observable quantities through the cross product which is defined 
so that it compensates for the inversion of pseudo vectors in absolute space. 

Such situation is not quite satisfactory however. Although many quantities in 
classical physics are not directly observable (measurable) and hence their 'realism' can 
always be questioned by a radical empiricist, a natural trend is to put as much 'realism' 
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into each quantity as possible. In this paper we discuss an alternative convention for 
the cross product that allows to eliminate pseudo quantities from classical physics. 
Interestingly, it turns out to be more than just a matter of convention or convenience 
once the complex-coordinate transformations are concerned. 

2. Maxwell's equations 

Classical theory of electromagnetism is usually built upon the fundament of Maxwell's 
equations — a set of dynamic field equations which we write in standard three- vector 
notation as 

V*E=-B, V-B = 0, (1) 

V*H = b+j, V-D = p, (2) 

accompanied, in the simplest case, by the constitutive relations D = eE and 
B = /il^J. Dots over characters denote time derivatives; tildes over H and B 
show explicitly that these quantities are regarded as pseudo vectors (such notation 
follows J. A. Schouten [T]); curls (V x ) lfe = ^ k dj are also defined through the pseudo 
density (£ 4:,fe , k = 1,2, 3, equal to the Levi-Civita permutation symbol e* jfe in any 
coordinate system [31 p. 158]: 

e'3' k ' = A- l Ji'j{4"e jk , e 3k = e ijk . (3) 

The = sign means equality in a given (here, in an arbitrary right- or left-hand) 
coordinate system. With use of classic 'kernel- index' notation, we write Maxwell's 
equations ([T]), in a more general form, valid assumably in arbitrary curvilinear, 
nonorthogonal coordinates: 

<),!■:,.. = diVBi* = 0, (4) 

e^djHk = ±> l + j\ <W = p. (5) 

Here -B, and Hi are the covariant vector and pseudo vector; D 1 and 25* the 
contravariant vector density and pseudo density of weight +1, as reflected by 
Gothic kernels; f and p the contravariant vector and scalar densities. This form 
follows from the (here unquestioned) three-dimensional, generally covariant Maxwell's 
equations [H[3] 

2d[iEfl — —Bij, d[iB jk ] = 0, (6) 

d d S) ij =i) 4 + f, d{S l = P , (7) 

if one makes use of the dual equivalents Q5 Z = ^^ k Bjk and £?< = ^}ijk& k , as 
suggested e.g. by equations (2.35) and (2.37) in [3]. In ((6]) and {7|, the square 
brackets denote alternation, B^ = —Bji is the covariant bivector, and ffi — —fy 1 the 
contravariant bivector density. Assuming the constitutive relations are *B* = p^ Hj 
and D 1 = e l ^Ej, one concludes that the permittivity and permeability are tensor 
densities of weight +1, so that they are transformed as 

- |A|- X jfjf^, = a './:'./;'//'•'. (8) 

An evident alternative to (|4|) , ([5]) can be constructed by using in ((6]) and {?]) the duals 
W = \<£ l 3 k Bjk and Hi = \tijk& k instead, with the ordinary densities € y7c and 

% To end up with SI units, one puts e=ee and u =A 1 M with e= 10 7 /(47rc 2 ) and M= 47rl0 -7 . 

to r ' v ' 
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equal to the Levi-Civita symbol for proper transformations, and changing the sign for 
improper ones: 

e 'j'k' = \A\- 1 Ji'jfjjf& k , e jk = e ijk . (9) 

Here = stands for equality in an arbitrary but right-hand coordinate system. We see 
no reason to agree that the duals constructed with <£ yfc and tijk "are only meaningful 
for proper transformations" [51 p. 41]. An appropriate form of the generally covariant 
Maxwell's equations — compare to ([5]) — reads 

e jk d 3 E k = -SB*, <W = 0, (10) 

<£ ijk djH k = & + f , (9,2)* = p. (11) 

Given the constitutive relations are SB 1 = p^Hj and = e'^Ej, the transformation 
laws for the permittivity and permeability are exactly the same as given by ©. In 
standard three- vector form, (|10p . (jlip are reduced to the commonly known equations 
HJ), © except for the tildes: 

V x E = B, VB = 0, (12) 

VxH = D+j, VD = p. (13) 

Two differences from the standard formulation thus is that all vectors are 'polar', and 
the curl VxF = ^ k djF k . This seems perfectly sound were we to assign any 'rigidity' 
to the picture of magnetic 'field lines' for stationary media, instead of absolutizing the 
cross product definition by use of pseudo permutation objects. 



3. Equivalence in real space 

Most often in formulas that define observable (directly measurable) quantities in terms 
of H or B, these latter fields enter in some conjunction with the vectorial product 
operation — like V x H in Ampere's law, j x B in the expression for the Lorentz force 
density, or H x in the terms in constitutive relations accounting for induced optical 
activity. If we note that the neighboring tildes, if present, one over cross product and 
other over the H or B vector, do always annihilate, we conclude that the 'polar-axial' 
and 'pure polar' formulations of electrodynamics are largely equivalent. 

Now we also show the equivalence of the two formulations in classifying the 
eigenmodes of symmetric systems by symmetry arguments. In the following, consider 
reflection with respect to the a; 1 = surface, x 1 — —x 1 . The corresponding Jacobian 
is J\ — diag(— 1, 1, 1) = J?,, its determinant A = —1. It is well known that reflection 
symmetry of the system, e'(x') = e(x) and /i'(x') = /x(x), allows to replace the whole 
system by its half, with the PEC (the 'electric wall') or PMC (the 'magnetic wall') 
boundary standing in place of the reflection symmetry plane. 

We write ((TDJ), (JTTJ) in terms of covariant electric and magnetic vectors: 

e 3k d 3 E k = -fi ij Hj, di^Hj = 0, (14) 

e ]k d 3 H k = e ij Ej + )\ die ij Ej = p. (15) 
In the 'pure polar' formalism, both vectors are transformed under reflection as 
E i ,(x')=4E i (j(x i ) 

= {-E 1 (-x 1 ,x 2 ,x 3 ),E 2 (-x\x 2 ,x 3 ),E 3 (-x\x 2 : x 3 )}, (16) 
H il ( X ') = 4H i (jfx i ) 

= {-i/i(-x 1 ,x 2 ,x 3 ),i7 2 (-a; 1 ,a; 2 ,x 3 ), J ff3(-x 1 ,a; 2 ,a; 3 )}. (17) 
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At the same time, if material objects are invariant under the given transform, so 
that e l 'j"(x') = e y (x) and // J "(x') = ^(x), then Maxwell's equations (H} and ([TB]) 
remain unchanged except for the —1 multipliers acquired by the left-hand sides, owing 
to the transformation rule for € y ©• Thus the admissible new solutions (in the 
absence of free charges) are B^(x') = ±aEi{x.) and 7Ji/(x') = ^aiJi(x). Invariance of 
the physical system (material objects plus electromagnetic fields) to double reflection 
gives a 2 = 1. Equating (Til))) with Ei and (fTT|) with —Hi yields 

- E^-x^x 2 ^ 3 ) = E^x^x 2 ^ 3 ), -H^-x 1 ^ 2 ^ 3 ) = -H^x 1 ^ 2 ^ 3 ), (18a) 
E 2 (-x\x 2 ,x 3 ) = E 2 (x\x 2 ,x 3 ), H 2 (~x\x 2 ,x 3 ) = -H 2 (x\x 2 ,x 3 ), (186) 
^(-a; 1 ,^ 2 ,^ 3 ) = ^(a; 1 ,^ 2 ,^ 3 ), H 3 (-x\ x 2 , x 3 ) = -H 3 (x\ x 2 , x 3 ). (18c) 

These conditions define even modes with the 'numerical' PMC plane at x 1 = 0: 
Ei\x 1 =q = 0, H 2 \ x i =0 — 0, H 3 \ x i =0 = (the 'physical' perfect conductor planes 
are not so readily expressed in nonorthogonal coordinates). The conditions obtained 
in a similar way with a = — 1 define the numerical PEC boundary and odd modes. 
Thus, the 'pure polar' Maxwell's equations (fTU|) . ([Ti"j) lead to exactly the same mode 
classification in optical waveguides and resonators exhibiting reflection symmetry as 
does the standard 'polar-axial' formulation. 



4. Complex-coordinate scaling 

In computational electrodynamics, the dominating way for treating unbounded 
problems in the finite-difference time-domain (FDTD), finite-difference frequency- 
domain, or finite-element calculations is with the perfectly matched layer (PML) 
method. Introduced originally in the 'split-field' formulation by J.P. Berenger [4] and 
soon afterwards recast in the 'uniaxial' [5j|6] and stretched-coordinate [71 l8l l9l ITOl ITT] 
forms, the PML concept stands as "one of the most significant advances in the 
historical development of the FDTD method" [T2] . 

The method is closely related to complex-coordinate transforms in classical 
quantum theory of atomic resonances [13j — though this is rarely appreciated by 
the electromagnetics modelling community. In electromagnetics, the method is 
implemented with surprising ease. Without loss of generality, we consider domain 
truncation in one, x 1 = x, direction. For media described by diagonal e and \i matrices 
in the given coordinates, those matrices are modified in the PML regions according to 

/ l/s x \ 
e = eA, A = MA, A= s x (19) 

x \ s x J 

with the complex function 

Sx { x '^) = ax') + \^1 (20) 

where x' G [—6, b] is the bounded computation-space coordinate, C( x ') an( i a ( x ') 
are real- valued functions equal to 1 and, respectively, over [—a, a] C [—6,6], while 
cr{x') > for a < \x'\ < 6 in order to damp the exp[i(kx' — ujt)] oscillatory waves inside 
the PML regions. Odd-power frequency dependence is included in the imaginary part 
of ((2T))) . and hence in (Til?]) , for causality reasons [HI §123]. It is repeatedly claimed that 
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Figure 1. Mode index errors on Cartesian (left) and polar (right) grids. The 
superiority of the standard perfectly matched layers compatible with the 'pure 
polar' Maxwell's equations dlOl . lilt is evident. 



the structure of permittivity and permeability modified within the PMLs according 
to (|19[) can be recovered by complex coordinate scaling 

x = I s x (x', us) da/. (21) 
J o 

It is easy to see however, that only the 'pure polar' formulation (fTUf . (fTTj) leads to 
1(15)). In contrast, © and © yield 

/ kx|/4 o o \ 

A= Is, I (22) 
x V \s x \ J 

which is reduced to (fT9| for real- valued coordinate squeezing [15] , but essentially differs 
from (| 19(1 given the s x function complex. 

In order to compare the performance of ([22| and (fT9|) . we simulated the 
fundamental mode of a step-index fiber (n cor o/n c iad = 1.515/1.5 as in [El [16]) with 
the two-dimensional finite-difference frequency-domain method on Cartesian and polar 
grids. In figure[T] the fundamental mode index errors are plotted, associated with the 
standard ('pure polar' consistent) PMLs, the 'polar-axial' PMLs, and simple zero 
boundaries. The corresponding physical domains are [0...1.5i?] x [0...1.5i?] and 
[0 . . . 1.5R] x [0 . . . 7r/2], with resolution 40 x 40 and 40 x 20 pixels, respectively. 

5. Summary 

Standard ('polar-axial') and alternative ('pure polar') formulations of Maxwell's 
equations were presented in section [2] In most cases, as illustrated in section [3[ both 
approaches yield equivalent description of electromagnetics phenomena. In section [4] 
we pointed however that it is the 'pure polar' formalism that complies with the PML 
method used with much success in computational electrodynamics. An interesting 
question is, To which extent can the classic electric and magnetic fields — whose parity 
is a matter of controversy — be associated with the 'electric 2-'-pole' and 'magnetic 
2 J -pole' photon states of different parity in quantum electrodynamics? 
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